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We present results of reanalysis on proton decay in the minimal SU(5) SUSY 
GUT model. Unlike previous analyses, we take into account a Higgsino dressing 
diagram of dimension 5 operator with right-handed matter fields {RRRR opera- 
tor). It is shown that this diagram gives a significant contribution for p — » K~^Ut- 
Constraints on the colored Higgs mass Mc and the sfermion mass m ^ from Super- 
Kamiokande limit become considerably stronger than those in the previous analy- 
ses: Mc > 6.5 X 10l6 GeV for < ITeV. The minimal model with mj <, 2TeV 
is excluded if we require the validity of this model up to the Planck scale. 



1 Introduction 

Supersymmetric grand unified theory (SUSY GUT) ||] is strongly suggested 
by gauge couphng unification around Mx ^ 2 x 10^^ GeV In this theory, 
the hierarchy between the weak scale and the GUT scale AIx is protected 
against radiative corrections by supersymmetry. Also, this theory makes suc- 
cessful prediction for the charge quantization. Proton decay is one of the direct 
consequences of grand unification. The main decay mode p K^V [Q in the 
minimal SU(5) SUSY GUT model |^ has been searched for with underground 
experiments [^,0], and the previous results have already imposed severe con- 
straints on this model. Recently new results of the proton decay search at 
Super-Kamiokande have been reported |^ . The bound on the partial lifetime 
of the K^V mode is r(p — > K^V) > 5.5 x 10'^^ years (90% C.L.)[|, where three 
neutrinos are not distinguished. 

There are many analyses on the nucleon decay in the minimal SU(5) SUSY 
GUT model [p|,p^|-p^. In the previous analyses, it was considered that the con- 
tribution from the dimension 5 operator with left-handed matter fields {LLLL 
operator) was dominant for p — > K~^v In particular a Higgsino dressing 
diagram of the RRRR operator has been estimated to be small or neglected in 
these analyses. It has been concluded that the main decay mode is p K^v^, 
and the decay rate of this mode can be suppressed sufficiently by adjusting rel- 
ative phases between Yukawa couplings at colored Higgs interactions . 

"^Talk given by T. Nihei in International Symposium on Supersymmetry, Supergravity, and 
Superstring (SSS99), Seoul, Korea, June 23-27, 1999, based on the published work [y. 
''After we finished our analysis, the latest limit of the Super-Kamiokande t(p K'^u) > 
6.7 X 10'^'^ years (90 % C.L.) |b[ was reported. An adaptation to the updated experimental 
limit is straightforward (See Eq.{B|)). 
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In this talk, we present results of our analysis on the proton decay including 
the RRRR operator in the minimal SU(5) SUSY GUT model @]|. We calculate 
all the dressing diagrams (exchanging the charginos, the neutralinos and 
the gluino) of the LLLL and RRRR operators, taking account of various 
mixing effects among the SUSY particles, such as flavor mixing of quarks 
and squarks, left-right mixing of squarks and sleptons, and gaugino-Higgsino 
mixing of charginos and neutralinos. For this purpose we diagonalize mass 
matrices numerically to obtain the mixing factors at 'ino' vertices and the 
dimension 5 couplings. Comparing our calculation with the Super-Kamiokande 
limit, we derive constraints on the colored Higgs mass and the typical mass 
scale of squarks and sleptons. We find that these constraints are much stronger 
than those derived from the analysis neglecting the RRRR effect. 

2 Dimension 5 operators in the minimal SU(5) SUSY GUT model 

Nucleon decay in the minimal SU(5) SUSY GUT model is dominantly caused 
by dimension 5 operators []l0| , which are generated by the exchange of the 
colored Higgs multiplct. The dimension 5 operators relevant to the nucleon 
decay are described by the following superpotential: 



Here Q, If^ and E"^ are chiral superfields which contain a left-handed quark 
doublet, a charge conjugation of a right-handed up-type quark, and a charge 
conjugation of a right-handed charged lepton, respectively, and are embedded 
in the 10 representation of SU(5). The chiral superfields L and D'^ contain a 
left-handed lepton doublet and a charge conjugation of a right-handed down- 
type quark, respectively, and are embedded in the 5 representation. A mass 
of the colored Higgs superfields is denoted by Mc- The indices i,j,k,l = 
1,2,3 are generation labels. The first term in Eq. (|^) represents the LLLL 
operator which contains only left-handed SU(2) doublets. The second term 
in Eq. (|^) represents the RRRR operator which contains only right-handed 
SU(2) singlets. The coefficients Csl and C5U in Eq. (|l]) at the GUT scale are 
determined by Yukawa coupling matrices as 




(1) 




iYDh{V^PYuV)ki, 
{P*V*YDhiYuV)ki, 



(2) 



=See also Ref. Q. 
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where Yu and Yd are diagonalized Yukawa coupling matrices for 10 • 10 • 5h 
and 10 • 5 • 5h interactions, respectively. The unitary matrix V is the Cabibbo- 
Kobayashi-Maskawa (CKM) matrix at the GUT scale. The matrix P — 
diag{Pi , P2 , P3) is a diagonal unimodular phase matrix with \Pi\ = 1 and 
detP = 1. We parametrize P as 

P1/P3 = e*'^i-% P2/P3 ^ e"^^\ (3) 

The parameters (jji^ and (j)23 are relative phases between the Yukawa couplings 
at the colored Higgs interactions, and cannot be removed by field redefinitions 
| p6[ . The expressions for C^l and C^r in Eq. (j^) are written in the flavor basis 
where the Yukawa coupling matrix for the 10 • 5 • 5h interaction is diagonalized 
at the GUT scale. Numerical values of Yjj, Yd and V at the GUT scale are 
calculated from the quark masses and the CKM matrix at the weak scale using 
renormalization group equations (RGEs). 

In this model, soft SUSY breaking parameters at the Planck scale are de- 
scribed by mo, Mgx and Ax which denote universal scalar mass, universal 
gaugino mass, and universal coefficient of the trilinear scalar couplings, re- 
spectively. Low energy values of the soft breaking parameters are determined 
by solving the one-loop RGEs |jl^. The electroweak symmetry is broken ra- 
diatively due to the effect of a large Yukawa coupling of the top quark, and 
we require that the correct vacuum expectation values of the Higgs fields at 
the weak scale are reproduced. Thus we have all the values of the param- 
eters at the weak scale. The masses and the mixings are obtained by diag- 
onalizing the mass matrices numerically. We evaluate hadronic matrix ele- 
ments using the chiral Lagrangian method fl^ ]. The parameters ap and Pp 
defined by (0|e''^'=(d5jU^X|p) = apNL and (0|e'''"=(d|w|,X|p} = /3pNL (Nl 
is a wave function of the left-handed proton) are evaluated as 0.003 GeV^ < 
Pp < 0.03 GeV^ and ap = -f3p by various methods [^9|. We use the smallest 
value Pp = —ap = 0.003 GeV^ in our analysis to obtain conservative bounds. 
For the details of the methods of our analysis, see Ref. Q|. 

3 RRRR contribution to the proton decay 

The dimension 5 operators consist of two fermions and two bosons. Eliminating 
the two scalar bosons by gaugino or Higgsino exchange (dressing), we obtain 
the four-fermion interactions which cause the nucleon decay Q,^]. In the 
one-loop calculations of the dressing diagrams, we include all the dressing 
diagrams exchanging the charginos, the neutralinos and the gluino of the LLLL 
and RRRR dimension 5 operators. In addition to the contributions from the 
dimension 5 operators, we include the contributions from dimension 6 operators 
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Figure 1: Higgsino dressing diagram which gives a dominant contribution to the p — > K^Vt 
mode. The circle represents the RRRR dimension 5 operator. We also have a similar diagram 
for {uRSR){dLUrL)- 



mediated by the heavy gauge boson and the colored Higgs boson. Though 
the effects of the dimension 6 operators 1/M|-) are neghgibly smaU for 
p — > K~^v, these could be important for other decay modes such as p — > 7r°e^. 
The major contribution of the LLLL operator comes from an ordinary diagram 
with wino dressing. The major contribution of the RRRR operator arises 
from a Higgsino dressing diagram depicted in Fig. 11^. The circle in this figure 
represents the complex conjugation of C^"^' in Eq. (H) with i ^ j ~ \ and k — I 
— 3. This diagram contains the Yukawa couplings of the top quark and the tau 
lepton. The importance of this diagram has been pointed out in Ref. jl5| in the 
context of a SUSY SO(IO) GUT model. The contribution of this diagram has 
been estimated to be negligible or simply ignored in previous analyses in the 
minimal SU(5) SUSY GUT fl | l0| -|ll. In particular, the authors of Ref. ^ 
calculated the diagram in Fig.^. However the effect was estimated to be small, 
because the authors assumed a relatively light top quark {mt ~ 50 GeV) . We 
use an experimental value of the top quark mass, and show that this diagram 
indeed gives a significant contribution in the case of the minimal SU(5) SUSY 
GUT model also. 

Before we present the results of our numerical calculations, we give a rough 
estimation for the decay amplitudes for a qualitative understanding of the 
results. In the actual calculations, however, we make full numerical analyses 
including contributions from all the dressing diagrams as well as those from 
dimension 6 operators. We also take account of various effects such as mixings 
between the SUSY particles. 

Aside from the soft breaking parameter dependence arising from the loop 
calculations, relative magnitudes between various contributions can be roughly 
understood by the form of the dimension 5 operator in Eq. (H). Counting the 
CKM suppression factors and the Yukawa coupling factors, it is easily shown 
that the RRRR contribution to the four-fermion operators {uiidR){sLVTL) and 
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{uBSR){dLVTL) is dominated by a single Higgsino dressing diagram exchang- 
ing (the right-handed scalar top quark) and f/j (the right-handed scalar tau 
lepton). For K'^v^ and K'^Ve, the RRRR contribution is negligible, since it 
is impossible to get a large Yukawa coupling of the third generation without 
small CKM suppression factors in this case. The main LLLL contributions 
to {uLdL){sLViL) and {uLSL){dLi'iL) consist of two classes of wino dressing 
diagrams; they are cl exchange diagrams and II exchange diagrams jl^. Ne- 
glecting other subleading effects, we can write the amplitudes (the coefficients 
of the four-fermion operators) for p — > K^Vi as, 

Amp.(p -> K+Ve) - [P2A^{£l) + Pj.Ae{tL)]LLLL, 

Amp.(p ^ K+V^) - [P2Am(5l) + PzA^{iL)]LLLL, 

Amp.(p K^Vt) ^ [-P2^t(cl) + -P3^t(?l)]llll + [-Pl^r(?fl)] (4) 

where the subscript LLLL {RRRR) represents the contribution from the LLLL 
{RRRR) operator. The LLLL contributions for At can be written in a 
rough approximation as ^r(cL) ^ g'^YcYbV*^VcdVcsm^ / {Mcmj) and Ar{iL) 

^ 9'^YtYbV*^VtdVtsmy^^ /{Mcmj), where g is the weak SU(2) gauge coupling, 

and is a mass of the wino W . A typical mass scale of the squarks and the 
sleptons is denoted by mj. For Af^ and A^, we just replace Yi,V*fj in the ex- 
pressions for At by YsV*^ and Y^V*^, respectively. The RRRR contribution is 
also evaluated as v4^(tfl) ~ YdY^YrV^lVudVtsf^/ {Mcm^) , where fi is the super- 
symmetric Higgsino mass. The magnitude of ^ is determined from the radia- 
tive electroweak symmetry breaking condition, and satisfies > {m^^l in the 
present scenario. Relative magnitudes between these contributions are evalu- 
ated as follows. The magnitude of the cl contribution is comparable with that 
of the tL contribution for each generation mode: |v4i(ci)| ~ |v4i(ti)|. There- 
fore, cancellations between the LLLL contributions P2Ai{cL) and P^Ai{iL) 
can occur simultaneously for three modes p — *■ K^Vi {i — e, ^ and t) by ad- 
justing the relative phase (j)23 between P2 and P3 [0. The magnitudes of the 
LLLL contributions satisfy \P2Af,{cL) + PaA^ih)] > \P2Ar{cL) + P3Ar{h)\ > 
\P2Ae{cL) + P3Ae{iL)\ independent of 023- On the other hand, the magnitude 
of Ar{iji) is larger than those of Ai{cL) and Ai(fi), and the phase dependence 
of Pi^T-(ifl) is different from those of P2^i(cL) and P3Ai(fi). Note that Ai{cL) 
and Ai{i]^) are proportional to l/(sin/?cos /3) — tan (3+1/ tan/3, while AT{iji) 
is proportional to ^ (tan f3 + 1/ tan/3)^, where tan/3 is the ratio of the vacuum 
expectation values of the two Higgs fields. Hence the RRRR contribution is 
more enhanced than the LLLL contributions for a large tan /3 . 
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Figure 2: Decay rates r(p K'^Ui) {i = e, fj, and r) as functions of the phase ij>23 for tan/3 
= 2.5 The other phase (/.13 is fixed at 210°. The CKM phase is taken as Sn = 90°. We 
fix the soft SUSY breaking parameters as mo = 1 TeV, Mgx = 125 GeV and Ax = 0. The 
sign of the supersymmetric Higgsino mass fi is taken to be positive. The colored Higgs mass 
Mq and the heavy gauge boson mass My are assumed as Mq = My = 2 X 10^'' GeV. The 
horizontal lower line corresponds to the Super-Kamiokande limit t{p — > K~^u) > 5.5 X 10^'^ 
years, and the horizontal upper line corresponds to the Kamiokande limit r(p — > K^p) > 
1.0 X lO''^ years. 



4 Numerical results 

Now we present the results of our numerical calculations For the CKM 
matrix we fix the parameters as V^s = 0.2196, Vet — 0.0395, |T4&/Kb| = 0.08 
and (5i3 = 90° in the whole analysis, where Si^ is a complex phase in the CKM 
matrix in the convention of Ref. |2^]. The top quark mass is taken to be 175 
GeV ||2l[]. The colored Higgs mass Mc and the heavy gauge boson mass My 
are assumed as Mc = My = 2 x 10^^ GeV. We require constraint on 6 57 
branching ratio from CLEO p2| and bounds on SUSY particle masses obtained 
from direct searches at LEP ||2|], LEP II Q and Tevatron 

Let us focus on the main decay mode p K^V. We first discuss the effects 
of the phases 013 and ^23 parametrizing the matrix P in Eq. In Fig. ^ we 
present the dependence of the decay rates T{p ^ K~^Vi) on the phase 023- As 
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an illustration we fix the other phase 0i3 at 210°, and later we consider the 
whole parameter space of 0i3 and (/)23- The soft SUSY breaking parameters 
are also fixed as toq — 1 TeV, Mgx = 125 GeV and Ax = here. The sign 
of the Higgsino mass /i is taken to be positive. With these parameters, all the 
masses of the scalar fermions other than the lighter t are around 1 TeV, and 
the mass of the lighter t is about 400 GeV. The lighter chargino is wino-like 
with a mass about 100 GeV. This figure shows that there is no region for 
the total decay rate r(p ^ K^V) to be strongly suppressed, thus the whole 
region of 023 in Fig. is excluded by the Super-Kamiokande limit. The phase 
dependence of r(p K^Vr) is quite different from those oi T{p K^V^) 
and Tip K+Ve). Though r{p K+Vf,) and r{p K+V^) are highly 
suppressed around (/)23 ~ 160°, T(ji K^Vr) is not so in this region. There 
exists also the region 023 ^ 300° where r(p K^Vr) is reduced. In this 
region, however, T{p — > K'^V^) and T(j) — > K^Vf.) are not suppressed in turn. 
Note also that the K^Vr mode can give the largest contribution. 

This behavior can be understood as follows. For 77^ and Fg, the effect 
of the RRRR operator is negligible, and the cancellation between the LLLL 
contributions directly leads to the suppression of the decay rates. This can- 
cellation indeed occurs around 023 160° for both and Ve simultaneously 
in Fig. 1^. For Vr-, the situation is quite different. The similar cancellation 
between P2^t(cl) and P-^AT-itL) takes place around 023 ~ 160° for Vr also. 
However, the RRRR operator gives a significant contribution for Vt- There- 
fore, T{p K^Vr) is not suppressed by the cancellation between the LLLL 
contributions in the presence of the large RRRR operator effect. Notice that it 
is possible to reduce r(p — > K^Vr) by another cancellation between the LLLL 
contributions and the RRRR contribution. This reduction of r(p K^Vr) 
indeed appears around 023 ^ 300° in Fig. ^ The decay rate r{p — > /C+F^) 
is rather large in this region. The reason is that P2At{cl) and PsAritL) are 
constructive in this region in order to cooperate with each other to cancel the 
large RRRR contribution PiAr{tR), hence P2^^i(cl) and PsA^itL) are also 
constructive in this region. Thus we cannot reduce both T{p — > K~^Vt) and 
T{p K^V^) simultaneously. Consequently, there is no region for the total 
decay rate T{p K^V) to be strongly suppressed. In the previous analysis [|l3| 
the region 023 ~ 160° has been considered to be allowed by the Kamiokande 
limit t{p -> K+V) > 1.0 x 10^^ years (90% C.L.) [§. However the inclusion 
of the Higgsino dressing of the RRRR operator excludes this region. We also 
examined the whole region of 0i3 and 023 with the same values for the other 
parameters as in Fig. ||. We found that we cannot reduce both T{p — > K'^Vr) 
and T{p — > K^Vf^) simultaneously, even if we adjust the two phases 0i3 and 
023 anywhere. 
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Figure 3: Lower bound on the colored Higgs mass Mc as a function of the left-handed scalar 
up-quark mass mf^^ Q]. The soft breaking parameters mo, Mgx and Ax are scanned within 
the range of < mo < 3TeV, < Mgx < 1 TeV and -5 < Ax < 5, and tan/3 is fixed at 
2.5. Both signs of fj, are considered. The whole parameter region of the two phases 0i3 and 
023 is examined. The solid curve represents the bound derived from the Super-Kamiokande 
limit t(p — » K^v) > 5.5 X 10^^ years, and the dashed curve represents the corresponding 
result without the RRRR effect. The left-hand side of the vertical dotted line is excluded by 
other experimental constraints. The dash-dotted curve represents the bound derived from 
the Kamiokande limit on the neutron partial lifetime T(n — > K'^'U) > 0.86 X 10^^ years. 



Next we consider the case where we vary the parameters we have fixed 
so far. The relevant parameters are the colored Higgs mass Mc, the soft 
SUSY breaking parameters and tan/3. The partial lifetime r(p K^V) is 
proportional to AIq in a very good approximation, since this mode is dominated 
by the dimension 5 operators. Using this fact and the calculated value of 
t{p K^v) for the fixed Mc — 2x 10^® GeV, we can obtain the lower bound 
on Mc from the experimental lower limit on t(jp K^V). In Fig. we 
present the lower bound on Mc obtained from the Super-Kamiokande limit as 
a function of the left-handed scalar up-quark mass m^^ . Masses of the squarks 
other than the lighter t are almost degenerate with m^^ . The soft breaking 
parameters mo, Mgx and Ax are scanned within the range of < mo < 3 TeV, 
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< Mgx < ITeV and -5 < Ax < 5, and tan/? is fixed at 2.5. Both signs 
of fi are considered. The whole parameter region of the two phases 0i3 and 
(^23 is examined. The sohd curve in this figure represents the resuh with 
all the LLLL and RRRR contributions. It is shown that the lower bound 
on Mc decreases like ~ l/m^i^ as rriuj^ increases. This indicates that the 
RRRR effect is indeed relevant, since the decay amplitude from the RRRR 
operator is roughly proportional to ^/ {Mcroj) ~ l/(Mcmj), where we use 
the fact that the magnitude of /.t is determined from the radiative electroweak 
symmetry breaking condition and scales as ~ m The dashed curve in Fig. |^ 
represents the result in the case where we ignore the RRRR effect. In this case 
the lower bound on Mc behaves as ~ 1/ , since the LLLL contribution is 
proportional to m^/ (Mcto|) . 

The solid curve in Fig. |^ indicates that the colored Higgs mass Mc must 
be larger than 6.5 x 10^^ GeV for tan/? — 2.5 when the typical sfermion mass 
is less than 1 TeV. On the other hand, there is a theoretical upper bound of 
the colored Higgs mass Mc ^ 4 x 10^^ GeV from an analysis of RGEs when we 
require the validity of the minimal SU(5) SUSY GUT model up to the Planck 
scale 0. Then it follows from Fig. | that the minimal SU(5) SUSY GUT model 
with the sfermion masses less than 2 TeV is excluded for tan/3 = 2.5. The 
RRRR effect plays an essential role here, since the lower bound on to j would 
be 600 GeV if the RRRR effect were ignored. We also find that the Kamiokande 
limit on the neutron partial lifetime r(n K'^v) > 0.86 x 10"^^ years (90% 
C.L.) already gives a comparable bound with that derived here from the 
Super-Kamiokande limit on r(p — > K'^V), as shown by the dash-dotted curve 
in Fig. a 

FigTQ shows the tan (3 dependence of the lower bound on the colored Higgs 
mass Mc obtained from the Super-Kamiokande limit. Here we vary toq, Mgx 
and Ax as in Fig. ^. The phases (/)i3 and 023 are fixed as 013 = 210° and 023 
= 150°. The result does not change much even if we take other values of 0i3 
and 023- The region below the solid curve is excluded if m^i^^ is less than 1 TeV. 
The lower bound reduces to the dashed curve if we allow m^i^^ up to 3 TeV. 
It is shown that the lower bound on ALc behaves as ~ tan^ /3 in a large tan f3 
region, as expected from the fact that the amplitude of p — > K'^Vt from the 
RRRR operator is roughly proportional to tan^ (3/Mc- On the other hand the 
LLLL contribution is proportional to ^ tan /?/Mc, as shown by the dotted 
curve in Fig. ^. Thus the RRRR operator is dominant for large tan/? Jist . 
Note that the lower bound on Mc has the minimum at tan/? « 2.5. Thus 

Also it has been pointed out that there exists an upper bound on AIq given by Mc 
< 2.5 x 10^^ GeV (90% C.LJif we require the gauge coupling unification in the minimal 
contents of GUT superfields lisl . 
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Figure 4: The lower bound on the colored Higgs mass Mq obtained from the Super- 
Kamiokande limit as a function of tan/3 llj. The phase matrix P is fixed by <f>i3 = 210° and 
023 = 150° . The region below the solid curve is excluded if the left-handed scalar up-quark 
mass m^^ is less than 1 TeV. The lower bound reduces to the dashed curve if we allow rrtflj^ 
up to 3 TeV. The result in the case where we ignore the RRRR effect is shown by the dotted 
curve for ma, < ITeV. 



we can conclude that for other value of tan /3 the constraints on Ale Sind m j 
become severer than those shown in Fig. ^. In particular the lower bound on 
mj becomes larger than ^ 2 TeV for tan/? 2.5. 

The constraints obtained from the figures can be expressed as follows: 

Mc \ > /^T'='=P(p^ii:+77)y ( f3p 



6.5 X IQis GeV J ^ \5.5x lO^^years J \ 0.003 GeV^ 

for tan/3 « 2.5 



5.0 X 1017 GeV J ^ \ lO^^years J V 0.003 GeV^ 

for tan/3 ^ 5, 
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where t''^^{p — > K'^v) is an experimental lower limit for the partial lifetime of 
the decay mode p K^V. 

5 Conclusions 

We have reanalyzed the proton decay including the RRRR dimension 5 op- 
erator in the minimal SU(5) SUSY GUT model. We have shown that the 
Higgsino dressing diagram of the RRRR operator gives a dominant c:ontribii- 
tion for p — > -ftT+F^, and the decay rate of this mode can be comparable with 
that ofp ^ K^V^. We have found that we cannot reduce both the decay rate 
of p ^ K^Vt and that oi p —> K'^V^ simultaneously by adjusting the rela- 
tive phases ^13 and (/>23 between the Yukawa couplings at the colored Higgs 
interactions. We have obtained the bounds on the colored Higgs mass Mc and 
the typical sfermion mass rrij from the new limit on t{p — > K'^V) given by the 
Super-Kamiokande: Mc > 6.5 x IQi^ GeV for < 1 TeV. The minimal SU(5) 
SUSY GUT model with m j ^ 2 TeV is excluded if we require the validity of 
this model up to the Planck scale. 
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